clear 

load run1a_baseline_sr.mat

gamma = parms.gamma;
chi = parms.chi;

zetafun = @(N,phi) ((N + phi.*(1-N)).^(gamma - 1)).*...
    ((N + (1-N).*(phi.^(1-1/chi))).^((1-gamma)/(1-1/chi)));

Ngrid = linspace(0.9,1,101)';

FontSize = 22;
LineWidth = 2;

figure
hold on
plot(Ngrid, zetafun(Ngrid, 1),'-.b','LineWidth',LineWidth)
plot(Ngrid, zetafun(Ngrid, 0.6),'-r','LineWidth',LineWidth)
legend({['$\Phi_t=1$'],['$\Phi_t=0.6$']}, ...
    'Location', 'NorthEast', 'interpreter', 'latex', 'FontSize', FontSize)
legend boxoff

set(gca,'xlim',[0.9 1])
set(gca,'xtick',[0.9:0.02:1])
set(gca,'xticklabel',[0.9:0.02:1])
box on

% set(gca,'ylim',[1 1.1])

ax = gca;
ax.XAxis.FontSize = 16;
ax.YAxis.FontSize = 16;

xlabel('Employment, $N_t$', 'interpreter', 'latex', 'FontSize', FontSize)
ylabel('$\zeta_t$', 'interpreter', 'latex', 'FontSize', FontSize)